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, We explore the role of an intermediate state (phase) in homogeneous nucleation phenomenon 

by examining the decay process through a doubly-humped potential barrier. As a generic model 
we use the fourth- and sixth-order Landau potentials and analyze the Fokker-Planck equation for 
the one-dimensional thermal diffusion in the system characterized by a triple-well potential. In 
\^ . the low temperature case we apply the WKB method to the decay process and obtain the decay 

' rate which is accurate for a wide range of depth and curvature of the middle well. In the case 

of a deep middle well, it reduces to a doubly-humped-barrier counterpart of the Kramers escape 
rate: the barrier height and the curvature of an initial well in the Kramers rate are replaced by the 
' arithmetic mean of higher (or outer) and lower (or inner) partial barriers and the geometric mean of 

curvatures of the initial and intermediate wells, respectively. It seems to be a universal formula. 
In the case of a shallow-enough middle well, Kramers escape rate is alternatively evaluated within 
the standard framework of the mean-first-passage time problem, which certainly supports the WKB 
result. The criteria whether or not the existence of an intermediate state can enhance the decay 
, rate are revealed. 

a 

g 

o . 

The decay process for a metastable system received considerable attention for long time. By using the WKB analysis 
of the Fokker-Planck (FP) equation, van KampenQ investigated the role of fluctuations in thermal diffusion process 
| under a symmetric bistable system and found that the system's ergodicity enables a Brownian particle to escape over 
^ . the potential barrier. In asymptotic time regime, he also obtained the astonishingly accurate formula for the so-called 
Kramers escape rate, i.e., the rate at which a Brownian particle escapes from a potential well over the single- humped 
barrier. Later, Caroli et al.|2j explored the asymmetric bistable system using a path integral (instanton) formalism. 
The separation of the time scale required for the Kramers rate treatment was investigated by Kapral 3] by using the 
| projection operator approach. For a Langevin equation with spatially- and time-correlated colored noise in the case 
■ of weak damping, the staggered-ladder spectra of the corresponding FP equation are found£|. 

Recently, Nicolis and co-workers 1|| studied the decay in the system with a metastable intermediate state in the 
i context of the protein crystallization, and within the rate equation treatment they gave an important suggestion that 
the existence of the intermediate state can enhance the nucleation rate under some parameter range. We note that 
one-dimensional treatment of chemical reaction along the reaction coordinate in the free energy landscape is very 
I 1 useful to describe the above problem. 

In this paper, we investigate the role of an intermediate state in nucleation phenomenon based on the WKB 
analysis of FP equation for the one-dimensional order parameter, and present a doubly-humped-barrier counterpart 
of the Kramers rate0,Q. The decay rate is given as the 1st excited eigenvalue of the FP operator. As is well known, 
FP equation is transformed into the associated Schrodinger equation and in the case that a well-defined thermal 
equilibrium does exist, eigenvalues of FP equation are exactly the same as those of Schrodinger eauationpl IfjL ITU. 
Further, in the low-temperature case we can apply the WKB analysis to the system with a triple-well potential and 
5_J ' extend Kramers escape rate to the situation with a doubly-humped barrier. 

This paper is organized as follows. In Section 2 we reduce the FP equation into the Schrodinger equation and obtain 
the WKB quantization condition for a system with triple-well potential. Section 3 is concerned with the evaluation 
of the tunneling integral, leading to the extended Kramers formula for the doubly-humped potential barrier. In 
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FIG. 1: (a) Landau potential W(x) in the presence of intermediate state in the case of deep central cave with A = 2.5, /j, = 0.775; 
(b)Effective potential V(x) corresponding to Fig.l (a). 



Section 4 the condition for the enhancement of nucleation rate caused by the intermediate state is revealed. In 
Section 5 we numerically check the validity of our new formula for the nucleation rate. Appendix is devoted to the 
alternative evaluation of the nucleation rate in the case of a shallow middle well, within the standard framework of 
the mean-first-passage time problem. 



II. FOKKER-PLANCK EQUATION AND WKB ANALYSIS 

Relaxation process responsible for a homogeneous nucleation is described by FP equation as 



d „, , Id fdWtx) n d \ „, , 



Ot 



rj is the viscosity and 9 = ksT is the temperature. Hereafter we assume r\ = 1 for simplicity. In order to specify the 
problem, we consider the symmetric triple- well potential proposed by Nicolis et al. p|, 



W(x) = J dxx{x z - A)(V - n) (2) 

with A > n > 0, which has two potential minima at x = two maxima at x — ±^/JI and one local minimum at 

x = as shown in Fig. 1(a). This potential mimics the barrier accompanied with a cave on the top of the double- well 
potential, which represents the presence of the intermediate state. By the separation anzats P(x,t) = 4>(x)e^ Et , 
the FP equation is transformed into the associated Schrodinger equation, 

H<j>(x) = -9(/)" (x) + V(x)4>(x) = E(f>(x), (3) 

with the effective potential 

W n W" 

V(x) = — (4) 

w 46 2 w 

and 9 corresponding to 9 = in quantum mechanics with Planck constant h and mass m. Equation J3J has the 

self-evident lowest-eigenvalue Eq = with the eigenfunction (f>o(x) = e~^& and the remaining eigenvalues satisfy 
E n > Eq. The problem of the time evolution of the probability density P(x,t) initially distributed in the left well 
reduces to that of finding the lowest few eigenvalues of Eq. © . Eventually the decay rate T is given by the first-excited 
eigenvalue. 

As shown in Fig. 1(b), V(x) has also symmetric three wells and two barriers. We note: even if the central well is 
shallow in the original potential W(x) (see Fig. 2(a)), the effective potential V(x) is negative at the corresponding 

local minimum (V(0) = — W J * 1 < 0), and for the original triple- well potential the effective potential V(x) has three 
potential minima, all with negative values (Fig. 2(b)). Thus the lowest three eigenvalues of Eq.^jj including Eq = 
experience all the three potential minima and two maxima (see Figs. 1(b) and 2(b)). We are particularly interested 
in the situation where the lowest eigenvalues in each of three wells are almost degenerate without tunneling through 
barriers. Then, by applying the WKB method we shall remove the degeneracy by incorporating the tunneling effects. 
Let us denote the regions bordered by the classical turning points xq, ■ ■ ■ ,Xq as a, 6, c, d, and e (see Fig. 1(b)). The 



FIG. 2: (a)Landau potential W(x) in the presence of intermediate state in the case of shallow central cave with A = 2.5, \i = 0.3. 
(b)Effective potential V(x) corresponding to Fig. 2 (a). Note: although W(x) has a shallow central cave, the bottom of the 
middle well of V is negative unlike the case of double-well potential with a single-humped barrier. 



wave functions in each regions 4> a ,--- ,4> e are connected by the connection formula, 

— e ±*( s +i) (E>V)^— (e s ±V s ) (E<V) 
VP VP 2 

S=\j f pdx\, p=V2m\E~V(x)\ (5) 
11 Jyo 

where yo denotes a classical turning point. The wave function in the left-most region a, is written as 4> a — 
sin(^ J xo pdx + -j), since it does not contain a term which exponentially diverges in the region — oo < x < xq. By 
using the connection formula, we can successively connect the wave functions until <p e in the right-most region as 

4> e = —{(-2i cos S a cos S c e Mb+Md + - sin S a sin S c e- Mb+Md 

Vp 2 

+ icos > S' a sinS* c e Mi, - Md + ~ sin S a cos S c e- Mb - Md )e lS 'e~^ pdx e l % 
2 8 

+ (2i C o S S aC osS c e^ -isinS aSi nS c e-^ ^ cos S a sin S c e^ 

+ 1 smS a cos S c e- Mb - Md )e- lS *e^ K' pdx e -^}. (6) 
8 

S a and Mfj are the action and tunneling integrals in each well a = a, c, e and barrier = b,d, respectively. Hereafter, 
wherever necessary, we use the abbreviated notation a = a, c, e to indicate the wells a, c, and e, respectively. Since 
4> e in the right-most well should not contain a term, sin(i J^ 5 pdx — connecting to the term which exponentially 
diverges in the region x 5 < x < oo, we obtain the following relation as |9lll0| 

4 cot S a cot S c cot S e - cot S e e- 2M » - cot S a e~ 2Md - - cot S c e- 2{Mb+Md) = 0. (7) 

Due to the small transparency of each barrier, actions are approximately quantized as 5*0, ~ Tr(n a + i) and one sees 

cot S a w 7r fn a + - - So^J . (8) 

Since our interest lies in the low-lying energy levels, we consider the case n a = n c = n e = 0. In this case, for low-lying 
eigenvalue E, action is evaluated as S a (E) w (E — E a )^\ E=E + S(E a ) — n E ^J a and substituting © into 0, one 
obtains the equation in the form of cubic polynomials of eigenvalue E, 
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FIG. 3: Energy splitting and shifts in the case (i) of a deep middle well. The unperturbed three eigenvalues are almost 
degenerate. 

Herei? Q , V a = (E a — ^j^) andw Q = \{^\e=e ) _1 defined for each well a imply the "vacuum energy" , (approximate) 
potential minimum and the frequency, respectively. 

III. TUNNELING INTEGRALS AND DECAY RATE 

The precise solutions of Eq.© constitute the lowest three eigenvalues with the lowest one being the ground state 
energy E = 0. To reach the goal, the evaluation of tunneling integrals for each well is necessary. In the problem 
under consideration, the potential is assumed symmetric and thus V a — V e , Mb = Md, and uj a = Lo e - The tunneling 
integrals are perturbatively estimated with use of small parameter 9. The classical turning points gives the interval 
[SE3, X4] on which the tunnel integration is performed. As shown in Figs.l and 2, the interval contains only one local 
maximum at x = y of the original potential W . In this case the tunnel integral is calculated as follows: We divide 
the integral 



6M d = I \l^—-^dx 



(10) 



into three parts and perform the integration as 



W" 



BMd = I —( 1 - w^ + o{d))dx 

y+ A l-ow" w 12 ( w 12 x 2 
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y _ A v 2 v_ 40W" \ 29W"J "~ J y+A '2 
^M±^M + W(y) -eio g (-W»(y)V6) + 0(9). 



(1 



W' 2 



0{9 2 ))dx 



(11) 



We choose A = 



26 



so that the above-used Taylor expansions, yl + x = 1 



0(x 2 ) are possible. 



Equation @ is solved as follows: When the barrier heights are large enough, the terms including the exponentially- 
small tunneling factors are negligible on the left-hand side of 0, giving a set of the lowest eigenvalues at each of 



separate wells, E a , e = + ^ a an( l ^ c = + When the effect of the small transparency due to tunneling will 
be incorporated, these three levels show a splitting and shifts (see Figs. [31 and 0}. 

In such a case one can as yet suppress the fourth term on the left-hand side of J5J|, which is a product of two 
tunneling terms and is almost vanishing, and one obtains the exact eigenvalues: 

_hw a 

Ea — + Va, 



E+ = 



huj a 1 t/ _[_ huj c 

n i~ V a ~T n 



V, 



4V((^+K)-(^ + W + ^e-- d 



(12) 



where E- stands for Eq = 0. The decay rate is given as the difference of two low- lying eigenvalues, r = E a — E- 
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FIG. 4: Energy splitting and shifts in the case (ii) of a shallow middle well. The unperturbed eigenvalue for a central well is 
separated from the doubly-degenerate ones in the left-most and right-most wells. 

Depending on the depth of the middle well of V(x), there are two typical cases, 
caS e (a): |(%* + V a ) - + V a )\ « ^f^ e -M d . 
case (b): |(%- + K) - + K)| » ^"vTO e -^, 

The cases (a) and (b) corresponds to deep and shallow middle wells, respectively. The energy splitting and shifts 
in each of these cases are illustrated in Figs. Eland 01 Our principal interest lies in the case (a) that the unperturbed 
three levels without tunneling terms are almost degenerate. In this case the first term in the square root in Ea. lll2|) 
is negligible and the decay rate is dominated by the tunneling term as 



r _ r _ p / 4h 2 uj a u c e\W"(y)\ ^^ Arcslhh(1) „-W(y)-w { * 3) +w (v) -w M 

a V - 2 \W'(x 4 )W'(x 3 )\ [ ' 



with the turning points x 3 = y -fjj, x 4 = VA — yjJT^JT) an< ^ ^hc maximum point y = ^JJl in the tunneling integration. 
Using the harmonic approximation hiu a c _ e sa 29V " {x™ 1 ^) for the curvature of each well of V(x), (|13|) is rewritten 



in the lowest order of 6 as 

. ^2- 



-e 2 -^y^W''(0)W''(>/A)|W''( v /M)|e- (AUl+At/a)/29 (14) 



where AU\ = W(y/jl) — W(y/X) and AU2 = W(y/jl) — W(0) are the height of potential barriers at the top measured 
from the bottom of the left-most well and from the bottom of the middle cave, respectively. W"(0), W"(y/Ji) and 
W"(y/\) are the curvatures at the cave, at the barrier top and at the the bottom of the left-most well, respectively. 
These quantities are illustrated in FigEl In deriving l|14(l from (|13|l we used the expressions valid in the lowest order 

of 6, y/26V"(x aie ) = W"(VX), y / 28V"(x c ) = W"(0), W'(x 3 ) = y/26W"(0) and W'(x A ) = y / 2W(VA). 

The result in Ea. (|14|l is a doubly- humped-barrier counterpart of the Kramers escape rate for a single barrier, 
since the potential-barrier height and the curvatures of a well in Kramers rate are replaced by the arithmetic mean 
of higher(or outer) and lower(or inner) partial barriers (AU\ and AU2) and the geometric mean of curvatures of 
the initial and intermediate wells (W"(y/X) and W"(0)), respectively (see FigEJ. The result in Ea. 1)14(1 cannot be 
obtained within the standard framework of the mean-first-passage time problem. 

On the other hand, in the case (b), one can Taylor-expand the formula in Ea. ll2|l in the small tunneling term, 
finding E_ = E„ — j- 1 h e~ 2Md . In this case the decay rate is given by 

to a 2 |(™+K)-( £ f £ +V r c )| w 2 to J 



(V2 - l) 2 ^_ 2 ^ 2 y] W" (0)W" ( VA) I W"(y/p) I ^ _ £ _ L _^ 
^ + V a )-(^ + V c )\ 



E a -E_= ^ v " ; e 4 ' 2 ^ v r ' y ~' yvr ~" e-^F^. (15) 



Although in the case of a shallow middle well the difference of harmonic energy levels can not be written in terms of 
original potential, the activation energy, AU\ + AU2 with AU2 ~ is nearly equal to the barrier height at the barrier 
top measured from the bottom of the left-most well and the conventional Kramers escape rate is recovered. As is 
proved in Appendix, the result in Eq. i|15[l can be confirmed by alternatively solving the corresponding rate within the 
standard framework of the mean-first-passage time problem. 
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FIG. 5: The potential barriers AUi and AU2 and curvatures in the original potential W(x). 



IV. ENHANCEMENT OF NUCLEATION RATE 

We proceed to investigate the possibility of enhancement of nucleation process due to the existence of an intermediate 
state on the basis of the low-temperature formula for the doubly-humped barrier version of Kramers rate. 

In the case of a symmetric double- well potential in Fig|S| W(x) = \ — \vx 2 : the decay rate is calculated [l| in the 
same way as in the previous Section and is given by 



F' = — - \/W"(y/v)\W"(0)\e 



(16) 



with a potential barrier AU — W(0) — W(y/U) (see Fig|SJ). We note that a faster decay due to the existence of an 
intermediate state is guaranteed under the following conditions: 

(i) the Boltzmann factor for a doubly-humped barriers is larger than that for a single barrier; 

(ii) Provided that two kind of Boltzmann factors are identical, the pre-exponential factor for a doubly-humped 
barriers is larger than that for a single barrier. 

The condition (i) implies that the mean barrier height Au ^+ Au ^ for triple wells is less than the single-barrier height 
for the double-well AU. To be explicit, the condition (i) is given by 



^ >V /G(A^), G(\,n) = 



A 3 - 3A 2 /i + 6A^ 2 - 2fi 2 



(17) 



On the other hand, in the case (ii) the quasi-equality Au ^+ Au ^ ~ A£7 leads to 



v « y/G(^). (18) 

Namely, v is almost a function of A and /1. For the pre-exponential factors to satisfy T > T', we have 

A/it(A-/x)i > 2- 3 /V. (19) 

Equations H18fl and (|19fl satisfy the condition (ii). 

In the same way, enhancement of the nucleation due to the intermediate state is also expected in the case of the 



sixth-order double- well Landau potential, W(x) 



fx . For instance, the condition (i) is now given by 



k > (3G(A,/x))s. 



(20) 



V. NUMERICAL EVALUATION OF DECAY RATE 



To verify the results in the previous Sections, we shall investigate the decay rate numerically by calculating the 
second-lowest eigenvalue of the associated Schrodinger operator H = ^- + V(x) in Eq.@ 0,0]. By using the 
eigenvectors of the harmonic oscillator to geth er with the annihilation-creation operator representation of p and x, 
(n\H\m) is diagonalized straightforwardly |l4j . For the system with a triple- well potential, the decay rate is obtained 
numerically in the parameter range 9 — 0.1, 2 < A < 3, 0.3 < fi < 1, which is compared with the formula in Eq. i|14|) 
in the case of a deep middle well (Fig. 7). In Fig. 7(a), there certainly exists a coincidence region where the decay rates 



FIG. 6: Schematic symmetric double-well potential W(x) with (broken line) and without (solid line) a cave in the central 
barrier. AU implies a barrier height in the case of the double-well potential with no cave in the barrier. 
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FIG. 7: Comparison of the decay rate V between the analytic formula in Eg. 1141 and numerically obtained values (dark circles) 
in the parameter range 2 < A < 3, 0.3 < fi < 1 at = 0.1. (a) Coincidence region where the ratio of the numerical and 
analytical decay rates falls between 0.819 and 1.221. Along the coincidence line (solid-line), the symmetric level splitting is 
observed as depicted in Fig.lJ^J. (b)Comparison between numerical (dark circles)and analytical (solid-line) decay rates along 
the coincidence line in Fig. 7(a). 

r in both the numerical and analytical results (Ea. (|14Jl ) coincides (see Fig. 7(a)). The best agreement can be seen 
along the coincidence line. In the vicinity of the coincidence region, the curvatures of the left-most and the middle 
wells are almost the same in the original Landau potential W(x), and the harmonic levels for each of three wells in the 
effective potential V(x) are almost degenerate, satisfying the condition for the case (a) below Ea. (|12() . Thus one can 
see the doubly-humped barrier counterpart of Kramers formula to be justified near the coincidence line (see Fig. 7(b)) 
in the present context. 

In the case of a shallow central well lying far from this coincidence region, the degeneracy condition is no longer 
satisfied. At A = 3, /i = 0.3, for instance, the middle cave is astonishingly shallow. Here, instead of Ea. H14|) . the formula 
in Ea. (|15[l explains the numerical Boltzmann factor. In Appendix, we alternatively derive the doubly- humped barrier 
counterpart of Kramers rate as the inverse of the first passage time. Due to the saddle point approximation, however, 
this method is valid only in the case of a shallow central well. 



With use of the WKB method, we have analyzed the homogeneous nucleation phenomenon in systems with an 
intermediate state, and obtained the decay rate for the thermal diffusion over a doubly-humped barrier. The analytic 
result is applicable for a wide range of depth and curvature of the intermediate middle well. In the case of a deep middle 
well, in particular, the decay rate becomes a doubly-humped-barrier version of the Kramers escape rate: the barrier 
height and curvature of the initial well in the Kramers rate are now replaced by the arithmetic mean of higher (or 
outer) and lower(or inner) partial barriers and the geometric mean of curvatures of the initial and intermediate wells, 
respectively. We have confirmed the presence of the coincidence region (almost a narrow strip) in the parameter space 
for the Landau potential where this universal formula holds well, and also revealed the condition for the intermediate 
state to enhance the nucleation rate. In the case of a shallow middle cave far from the coincidence parameter region, 
however, we find a less-essential modification of the Kramers rate, which is also verified by an alternative study based 
on the standard framework of the mean-first-passage time problem . 



VI. SUMMARY AND DISCUSSIONS 
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FIG. 8: Schematic doubly- humped potential with a shallow cave, W(x). 
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APPENDIX A: METHOD OF THE MEAN-FIRST-PASSAGE TIME: THE CASE OF A SHALLOW 

MIDDLE CAVE 



For a triple-well potential W with doubly-humped barrier, the counterpart of Kramers escape rate is straightfor- 
wardly derived with the aid of the saddle point approximation for the mean- first-passage time|l6|. The left-most 
well is approximated as a harmonic potential and the middle well is given by the fourth-order symmetric potential, 
!|x 2 — jx 4 . The mean first passage time t(x) from an initial point at x is given as 



V W(») 



t(x) = rj / dye / dze « , (Al) 

J X J — oo 

where the boundary x a is larger than two potential maxima. We further assume: a Brownian particle can not come 
back once it leaves the domain [— oo, x a ]. At low temperature, the saddle-point approximation for the integral over z 
is performed as 



W(z) 



DC 



w(z) w min / 2irQ 



dze o w / dze ° =e » J— — . (A2) 

1 V ^min 



Here W m i n and Lo m in arc the potential minimum and the curvature at x m i n in the left-most well. 
The integration over y is performed as 

dye o w / dye 8 



= 2t b e e ^-iW+^Wj (A3) 

with W cave the potential minimum in the cave . The decay rate T is given as the inverse of r as 

I -^)+h(^)) ■ (A4) 



7ry2Wf 
of the as 

out to be the usual Kramers escape rate. 



In the low temperature, with use of the asymptotic form for the modified Bessel function Ik(z) ~ 7/§^j) Ea. ljA4|l turns 



9 



where AW = W cave — W m i n + stands for the net barrier height between the minimum of the left-most well and the 
potential top of the doubly- humped barrier and y/a corresponds to the curvature at the top, which verifies Ea. Q15(l . 
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